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Abstract 

We propose a double-well configuration for optical trapping of ultracold two-species Fermi-Bose 
atomic mixtures. Two signatures of macroscopic quantum coherence attributable to a superfluid 
phase transition for the Fermi gas are analyzed. The first signature is based upon tunneling of 
Fermi pairs when the power of the deconfining laser beam is significantly reduced. The second 
relies on the observation of interference fringes in a regime where the fermions are trapped in two 
sharply separated minima of the potential. Both signatures rely on small decoherence times for 
the Fermi samples, which should be possible by reaching low temperatures using a Bose gas as 
a refrigerator, and a bichromatic optical dipole trap for confinement, with optimal heat capacity 
matching between the two species. 

PACS numbers: 05.30.Fk, 32.80.Pj, 67.60.-g, 67.57.-z 
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I. INTRODUCTION 



Degenerate Fermi gases are rather ubiquitous in nature at both the microscopic and 
macroscopic level, from nuclear matter to neutron stars. Studies of both their non- interacting 
and interacting features allow for the understanding of a wealth of physical phenomena occur- 
ring in the mesoscopic realm, in particular superconductivity. More recently, the possibility 
to cool dilute samples of Fermi gases below the microkelvin range has opened up a novel 
route to identify some of the fundamental features underlying interacting many-body Fermi 
systems [l | . While Pauli blocking [3] and Fermi pressure 3, m have been already evidenced, 
:'ocus on interacting properties has recently led to interesting effects in Fermi-Bose mixtures 
ij], and in two-component Fermi gases p. In particular, in the latter case evidence has 
been reported for anisotropic free expansions of a Fermi cloud when this is brought to a 
deep degenerate regime. The data have been interpreted in terms of a superffuid state of 
the Fermi gas, as predicted in ?| . However, alternative interpretations in terms of hydro- 
dynamic behaviour of a high density Fermi gas are also plausible j^j]. More recently, various 
rave used resonant superfluidity P. lid llll . ll^ to explore the BEC-BCS crossover 



groups 
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14, 



15| . with various claims for the formation of bound states of many body character, 



as expected for instance by BCS-like couplings, based on the dynamics of formation of Fermi 
pairs and on collective properties This reminds of the previous situation of degenerate 
Bose gases, when various indirect evidences were collected for the existence of a macroscopic 
quantum state by studying collective properties. The final evidence was only achieved by 
explicitely showing quantum coherence [17[ and, thereafter, macroscopic quantum transport 
phenomena like superfluidity [3] and quantized vortices [3, Q . Analogously, we do expect 
the coherence of the macroscopic wave function associated to a Cooper-paired state of Fermi 
atoms to play an important role to assess its superfluid nature. In this paper, we discuss a 
configuration for an optical dipole trap that could allow for quantitative studies of quantum 
coherence in an ultracold Fermi gas. In Section II we describe a geometry for an optical 
dipole trap which creates a bistable potential for both the Fermi species and the Bose species 
necessary to sympathetically cool the Fermi gas. In Section III we discuss possible signa- 
tures for macroscopic quantum coherence through tunneling phenomena in a regime where 
the laser intensity of the blue-detuned beam is kept low. In Section IV we describe inter- 
ference experiments which should be able to disentangle between a BCS or a BEC regime 
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for the degenerate Fermi gas by observing the dynamics of the fringe visibility during the 
free expansion of the clouds. Macroscopic coherence in itself does not rely on the Fermi gas 
being in an effective BEC or a BCS state, as correlated Fermi pairs, either in a molecular 
state (BEC limit) or a many-body state (BCS limit) always behave as quantum coherent 
systems j^lj. However, for the interference fringe experiment and BCS-paired fermions, a 
sudden loss of fringe visibility is expected for large times of flight, while such a loss is not 
expected in the case of fermions coupled in a molecular state. Potential decoherence sources 
and some technical difficulties to be overcome are then discussed in the conclusive Section 
V. 



II. DOUBLE- WELL BICHROMATIC OPTICAL DIPOLE TRAPS 



The configuration we analyze relies on using an optical dipole trap made of focused red- 
detuned beams for trapping the atoms, and further blue-detuned beams for their selective 
deconfinement. Such a bichromatic optical dipole trap could allow to achieve a deep degen- 
erate regime for a Fermi gas when the latter is sympathetically cooled through a Bose gas 
undergoing evaporative cooling j^ . 

Let us consider an optical dipole trap consisting of a single red-detuned beam (optical 
source 1) propagating along the x-axis, and a single blue-detuned beam (optical source 2), 
also focused on the same spot, propagating along the orthogonal axis y. The resulting 
effective potential experienced by the atoms of species a (a = f for fermions and a = b for 
bosons) can be written as: 



U a (x,y,z) 
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where T" = l/(f2 Q — Qi) + l/(fl a + ^i) is a parameter related to the detuning between the 
atomic transition angular frequencies Q a = 27rc/A Q and the laser beam angular frequencies 
fij = 2-Kc/Xi (A a and Aj being the atomic transition and laser beam wavelengths, respec- 
tively), Pi and Wi are power and waist of the laser beams, Ri = nwf/Xi their Rayleigh 
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ranges, and P a at = hQ a T a /12irc 2 is the saturation intensity for the atomic transition. The 
potential (JIJ is well approximated by a second order Taylor expansion around the x axis 

U a ( Xl y, z) ~ U a (x) + l -d 2 y U a {x)y 2 + ^d 2 z U a (x)z 2 (2) 

where x = (x, 0, 0) and 



UJx) - 



4ttI 



sat 



w 2 



w 2 l + x 2 /R 2 

(3) 



+ exp (-2x 2 /w 2 2 ) 



The expression for U a (x) explicitly shows that there is a soft attractive potential on the 
Rayleigh range scale Ri and a sharp repulsion around the origin on the beam waist length- 
scale u>2- The net effect of the beams is to establish a double well potential along the x-axis 
with minima at ±x m = (±x m , 0, 0) having at the same time a strong quasi harmonic con- 
finement in the y — z plane. The transverse angular frequencies u ay (x) = (d 2 U a (x) /m a ) l l 2 
and u az (x) = (d 2 z U a (x) /m a ) l l 2 are one order of magnitude larger than the intra-well lon- 
gitudinal angular frequency u ax (x m ) = (d 2 U a (x m ) /ma) 1 ! 2 . For atomic gases with chemical 
potential fi a satisfying hu ax <C \i a hu ay , fuj az the present configuration thus realizes a 
quasi-lD trapped gas analogously to that obtained in highly elongated magnetic traps 
In Fig. Q]we show the potential energy for the Fermi and Bose components in the case 



2J| , by assuming a Nd: YAG 



of the 6 Li- 23 Na mixture, already brought to degenerate regime 
laser emitting at Ai = 1064 nm as red-detuned source, and its second harmonic as blue- 
detuned source. Details of the potential energy profiles are shown for the two species in 
Fig. |2j It is evident that the bosonic species experiences a double well potential with larger 
distance and higher potential barrier between the minima with respect to the fermionic one. 



Some comments are in order. The Fermi gas is always more deeply and strongly confined 
than the Bose gas, which is favorable for different reasons. Firstly, this allows for a continuous 
evaporative cooling of the Bose gas marginally affecting the Fermi gas. Secondly, the ratio 
between the average trapping frequency for fermions and bosons maintains a value greater 
than unity, with a negligible decrease with respect to the case of coaxial beams. This allows 
for efficient sympathetic cooling of a degenerate Fermi gas and a less degenerate Bose gas, 
through matching of the specific heats. Indeed the Bose gas, less degenerate, maintains 
its large classical specific heat at lower temperatures with respect to the more degenerate 
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FIG. 1: Bistable potential for optically trapped Fermi-Bose mixtures. The plots represent the 
equipotential surfaces above the minima of the potential f/™ m by an amount AU a =0.5, 2.5, 5, 10, 
20 nK (from inner to outer shells, respectively) for fermionic 6 Li (top) and bosonic 23 Na (bottom). 
We assume a laser power of Pi = 10 mW at Ai=1064 nm, Pi/P\ = 2.5 x 10~ 3 at A2=532 nm, and 
waists w\ = W2 = 10 /mi. The atomic transition wavelengths are Af = 671 nm and Ab = 589 nm. 

Fermi gas 2sj|. Thirdly, this reduces the peak density of the Bose gas at the center of 
the trap making less relevant both the boson-boson interaction, and most importantly the 
fermion-boson interaction. 

Finally, the fact that bosons experience a higher potential barrier and a larger separation 
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FIG. 2: Equipotential profiles in the x — y plane of U a (x,y,0) (left), for AU a =0.5, 2.5, 5, 10, 20, 
30, 40, 50, 60, 70 nK (from inner to outer shells, respectively), and potential energy along the x 
axis for the two species (right). The minima of the bistable potentials have been shifted to zero 
for the sake of comparison, their values being [7 f mm = —3.87 fiK and [/™ m = —3.31 fiK. All trap 
parameters as in Fig. ^ 

between the minima, strongly differentiates the dynamics of the two species. As explained in 
the following, the latter feature allows for two unambiguous signatures of the macroscopic 
coherence associated to a possible superfluid phase of the Fermi gas, based on tunneling 
oscillations and interference, respectively. 



III. TUNNELING PHENOMENA 



Macroscopic tunneling phenomena have been successfully explored in a Fermi liquid, 
namely 3 He j^J, and it is therefore natural to explore their counterpart in dilute Fermi 
gases. In our configuration, by using sufficiently low barriers the tunneling probability for 
fermions may become large enough to make the detection of the corresponding interwell 
oscillation feasible. Here there is a rich scenario due to the possibility of single-particle 
tunneling for both fermions in a degenerate but normal state and bosons in a thermal state 
(either quantum or thermally activated), and macroscopic quantum tunneling of fermions 
in a BCS state and of bosons in a BEC state. An oscillating tunneling current between the 
two wells of the trap is obtained by breaking the symmetry along the x axis by means of 
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a tilting potential, V(x) = b x, suddenly added to the trapping potential U a (x,y, z), e.g. 
using Zeeman shifts generated by a quadrupole magnetic field with symmetry axis along the 
x direction. 

The evaluation of the macroscopic tunneling current is particularly simple within the 
quasi one-dimensional approximation (J2J) of the potential. Suppose that the total number 
N a of trapped (bosons or bosonized fermions) atoms available for macroscopic tunneling 
is in the ground state of the potential Q and let ipa{x) be the unit-normalized associated 
wavefunction. If a tilting potential V(x) = b x with small bias parameter b is added at time 
t = 0, the evolution of the system wavefunction can be obtained as 

i/> a (x, t) = c^e lE " t/h <P1(x) + c%e iE % t/h ^{x), (4) 

where 4>i{x) and </>f (x) are the first two eigenf unctions, with eigenvalues and E% , re- 
spectively, of the tilted potential U a (x) + V(x) and the coefficients cf and cf are given by 

+00 

<j>?(x)i(; a (x)dx, i = 1,2. (5) 

■00 

In the absence of decoherence phenomena the evolution of the macroscopic wavefunction (j3J) 
gives rise to condensate fractions in the left and right wells, N E (t) and N^(t) oscillating in 
time with N%(t) + N^(t) = N a . The corresponding current is easily evaluated as 

~dt = a JtJ \Mx i t)\ 2 dx = A a an—^- } (6) 

with A a = 2c" cf Ci 2 N a AE a /h the amplitude of the tunneling current, depending upon the 
energy splitting AE a = E% — E", and 

r-O 
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cj>Ux)^(x)dx (7) 



is the overlap integral. 

To evaluate the macroscopic tunneling current for the realistic trap potential (J2J), the 
eigenvalues E" and E% and the corresponding eigenfunctions 0" and 02 m ust be found 
numerically. Due to the strongly differentiated tunneling dynamics for bosons and bosonized 
fermions, we have the hard numerical problem of finding exceedingly small energy splittings. 
The selective relaxation algorithm proposed in (27 1 allows to solve the problem at least in 
the parameter region of physical interest. 
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FIG. 3: Macroscopic coherence of the Fermi gas through tunneling experiments. Dependence of 
tunneling current amplitude per unit of available atoms A a /N a (solid line) and tunneling angular 
frequency AE a /h (dashed line) for 23 Na and 6 Li- 6 Li Cooper pairs versus bias strength b. For the 
sake of comparison, amplitude and frequency of bosons have been multiplied by 10 5 . The trap 
parameters are as in Fig. ^ 

In Fig. |3] we report the amplitude and the angular frequency of the currents for macro- 
scopic tunneling of 6 Li- 6 Li Cooper pairs and of the Bose condensed component of 23 Na versus 
the tilting parameter b. For both species there is an optimum value of the bias maximizing 
the observability of tunneling oscillations, with a maximum value A a ~ 4 s" 1 obtained for a 
bias value of the tilting potential b ~ 14 pK//im. Tunneling current and frequency for the 
Bose gas are smaller by five orders of magnitude with respect to the analogous quantities of 
the Cooper pair gas. This results from the exponential sensitivity of quantum tunneling to 
the different potentials experienced by the two species and to their different masses. 

Provided that the number of Cooper pairs is large enough, the modulation of the number 
of atoms in the two wells can be evidenced by using non-destructive phase-contrast imaging. 
The low intensity of the blue-detuned beam required to create a tunneling regime in the 
presence of a magnetic trap with limited confinement strength prevented observation of 
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tunneling phenomena in the experiment described in HH (see also |28( for a recently achieved 
bistable configuration). In the situation proposed here this issue is circumvented due to 
the possibility to change both the intensities of the beams while maintaining their ratio 
constant, at least to the extent that heating from residual Rayleigh scattering does not 
play a significant role. Current techniques allow for relative stabilization at the 1% level 
or below (for a general discussion of laser stabilization techniques see Q])' especially for 
frequency-doubled beams as in the proposed configuration. 

An important requirement for this proposed test is to maintain the decoherence rate low 
enough to minimize damping of the coherent oscillations which, according to Fig. 3, are 
expected to occur with periods of order 100 ms or longer. In the cooling strategy outlined 
in [2^ there is not a strict need to use enhancement of the elastic scattering length through 
Feshbach resonances 0, to reach a deep Fermi degenerate regime. This could circumvent 
the issue of decoherence sources due to enhancement of density, like those discussed in 
[30! for three-body collisions. Thus decoherence is mainly expected from the presence of 
Rayleigh heating or more technical sources like intensity and beam pointing fluctuations of 
the laser beams. Their effect is a temperature increase with a consequent increase of the 
thermal component and the single-particle tunneling current of bosons and fermions. 
It should be noted however that in an optical dipole trap, due to the smaller trapping 
volume and lower trap depth, we expect a suppressed thermal fraction. This is a further 
advantage in using an optical trap with respect to a magnetic trap for coherence experiments. 
Continuous evaporative cooling of the Bose species should also mitigate its effect. One can 
then envisage a cooling strategy where coherent oscillations persist for a much longer time 
at time-dependent amplitude and frequency, therefore originating a chirped signal. 

IV. INTERFERENCE PHENOMENA 

In a landmark experiment, the Ketterle group evidenced the macroscopic coherence of a 
pair of Bose condensates by looking at the interference fringes resulting after their release 
from a bistable potential 17[. The latter was obtained by the combination of a harmonic 
potential created by a magnetic trap, and a blue-detuned beam focused on the magnetic 
potential minimum with propagation orthogonal to the weaker confining axis. The average 
distance between the two condensates was controlled by changing the power of the blue- 
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detuned beam. 

A similar experiment can be repeated with an ultracold Fermi-Bose mixture in a bistable 
optical dipole trap. Above the temperature for the onset of a BCS-like phase transition 
of the Fermi gas one expects only interference fringes arising from the condensed fraction 
of the Bose gas. Below Tbcs we do expect also the emergence of an interference pattern 
coming from the bosonized fermions, a small fraction of the total number of Fermi atoms. 
The distance between the peaks of maximum signal in the interference pattern is given, for 
a free expansion, by l a = 2irht/m a d a . Here t is the time of flight, d a the initial distance 
between the centers of mass of the two clouds before the release from the trap, and m a 
either the mass of the Bose atoms or twice the mass of each fermion in the case of the 
BCS-bosonized component. Since the Fermi species in our example has both a significantly 
smaller mass and separation d a than the Bose species we do expect an easy discrimination of 
the interference pattern attributable to the former species. This is confirmed by looking at 
the distance between the interference peaks for 6 Li and 23 Na versus the P2/P1 power ratio, as 
depicted in Fig. 0] One can take advantage of this dependence of the interference patterns 
to discriminate any effect of the Bose component. The spacing of the interference fringes 
expected from the macroscopic wavefunction associated to the bosonized Fermi component 
maintains a value ~ 2 times larger than the corresponding one for the Bose s pec ies. Selective 
optical pumping tomography on the Fermi species cycling transition as in |l7[ can allow to 
enhance the corresponding interference signal. Considering the very dilute nature of the 
bosonized Fermi gas, as a consequence of the less stiff confinement of the Bose species, we 
do not expect a significant decrease of the fringe visibility due to the mean field effects as 
instead already observed for bosons j^. Due to the lower trapping frequencies when using 



a bichromatic optical trap as discussed in 



221, the mean-field effects due to the interaction 



between the Fermi and the Bose e also negligible. 

The interference experiment discussed above could also allow to distinguish if the Fermi 
gas is in a BCS or BEC regime. For BCS macroscopic states we do expect a loss of fringe 
contrast at expansion times much smaller than those reachable in the case of molecular 
BEC states. This behaviour is quantitatively shown by assuming an adiabatic expansion of 
the gas after the release from the trap, and a rough estimate of the critical temperatures, 
Tbcs an d Tbec as follows. Suppose that at time t — 0, when the trap potential is turned 
off, the gas made of Np fermions of mass m-p has density no- Let Rq be the radius of the 
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FIG. 4: Macroscopic coherence of the Fermi gas through interference experiments. Dependence 
of fringe spacing l a for 23 Na bosons and 6 Li- 6 Li Cooper pairs versus Pi/P\. We assume a time of 
flight t = 10 ms. 

equivalent sphere containing the gas at t — defined by |7ri?Qno = N F . At later times the 
cloud expands and, assuming an ideal behavior, the radius of the equivalent sphere is given 
by the law R(t) = Rq + Vpt. As a consequence, the density of the gas n(t) decreases from 
its initial value no and, in turn, the Fermi velocity Vp = Ufi^n) 1 ' 3 /m-p also decreases. The 
density of the gas is therefore determined by the following self-consistent equation 

*(*) = — ? " ( 8 ) 



which, for a generic value of t, must be solved numerically. For t^> to — m^Ro / h^ir 2 ^) 1 / 3 , 
the initial radius Rq can be neglected in the denominator of (jSJ) and we obtain the asymptotic 
time- varying density 



= (=*«»). (9) 

Once the density n(t) is known, we have an explicit estimate of the critical temperature for 
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FIG. 5: Temperature of the cloud (dashed line), and critical temperatures Tbcs> ^bec (continuous 
lines) as a function of time t in the case of a freely expanding 6 Li cloud with Np = 3 x 10 6 atoms. 
The initial density is no = 3.5 x 10 13 cm -3 and we use an elastic scattering length for fermions 
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l|. The dashed lines close to the critical temperature continuous lines are the 



asymptotic values for Tbcs> Tbec obtained by substituting n(t) with n^i) into (|TU|) and 
BCS transition jlj as a function of time 

71 



5 h\^n{m 
TbcsW = — — — — ex P 



2|a|(37rn(t))i 



3e 2mpkB 

where a is the elastic scattering length of the fermionic species. The critical temperature 
Tbec for fermions condensed via molecular states can be estimated from the critical tem- 
perature of an ideal gas of bosons with mass 2m F and density n(t)/2 



T (t)~ ^ ( H{t) 



(ii) 



where £ is the Riemann zeta function. 

In Fig. El we show the behavior of the critical temperatures Tbcs and Tbec> an d of the 
temperature of the cloud assuming an adiabatic expansion versus time in the case of 6 Li. 
While T B ec nas om y a n(t) 2 ^ 3 dependence, and therefore decreases for large times as t~ x , 
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Tbcs has a further exponential suppression factor and for large t decreases as t _1 exp (—\ft). 
By assuming an initial temperature T = 5 x 10~ 3 Tp, the disappearance of the fringes at 
times t ~ 10 ms when T/Tbcs > 1 would imply a BCS regime for the degenerate Fermi 
gas. On the other hand, the observation of interference effects should be possible up to 
longer times in the case of fermions coupled through molecular BEC, since the scaling of the 
temperature of the cloud and Tbec are similar. 



V. CONCLUSIONS 



We have discussed a quasi one- dimensional bistable configuration for optically trapped 
atoms. Two signatures have been discussed for evidencing macroscopic quantum coherence 
of a paired component of a Fermi gas, regardless of the BEC or BCS-like regime for the Fermi 
gas. Our proposal allows to identify a superfluid component both in a strongly-coupled 
regime based upon enhancement of scattering length (molecular BEC regime) or in a BCS- 
ike regime obtained by just cooling the sample at very low temperatures as suggested in 
25^ , with efficient heat capacity matching between the Fermi and the Bose species, without 
necessarily exploiting Feshbach resonances to obtain large critical temperatures for BCS 
pairing, although their use is certainly possible in our framework. The use of a Fermi-Bose 
mixture seems preferable since analogous quantum coherence experiments involving mixtures 
of two Zeeman levels of fermions are less easy to perform. Indeed, two distinguishable Fermi 
states with same mass will give rise to two independent interference patterns or tunneling 
currents just differing by the initial random phase, unless a locking mechanism is used. Also, 
dual evaporative cooling substantially decreases the number of atoms potentially available 
for Cooper pairing, then sensibly diminishing the expected signal. Moreover, the presence 
of the Bose gas until the last stage of the cooling is useful to quantitatively assess the 
temperature of the Fermi gas, and allows to study a variety of situations for which changes 
to the BCS-coupled Fermi pairs are expected when bosons mediates their interactions |34J. 

There have been many refined and ingenious proposals for the observation of the su- 
perfluid phase of an ultracold Fermi gas, ranging from the study of collective modes (^j], 

density profile ofthe Cooper-paired component light scattering 



moment of inertia 

210,0,0, 41 



42|, free expansion p|, Bloch oscillations 



44j . and Raman photoassociation in Bose- Fermi mixtures [4j|. All these proposals, includ 




internal Josephson effect 
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ing ours, will have to face the small number of fermions available in the superffuid state, 
and the subsequent small signal-to-noise ratio for any conceivable signature of the phase 
transition. In view of these experimental difficulties it is crucial to seek for redundancy 
of signatures with diverse techniques, hopefully all converging in individuating a common 
superffuid phase diagram. 
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